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ABSTRACT 

We show that the ratio of molecular to atomic gas in galaxies is determined by hydrostatic pressure 
and that the relation between the two is nearly linear. The pressure relation is shown to be good over 
three orders of magnitude for 14 galaxies including dwarfs, H I-rich, and H2-rich galaxies as well as 
the Milky Way. The sample spans a factor of five in mean metallicity. The rms scatter of individual 
points of the relation is only about a factor of two for all the galaxies, though some show much more 
scatter than others. Using these results, we propose a modified star formation prescription based on 
pressure determining the degree to which the ISM is molecular. The formulation is different in high 
and low pressure regimes defined by whether the gas is primarily atomic or primarily molecular. This 
formulation can be implemented in simulations and provides a more appropriate treatment of the 
outer regions of spiral galaxies and molecule-poor systems such as dwarf irregulars and damped Lya 
systems. 

Subject headings: GalaxiesdSM — ISM:evolution — ISM:molecules 



1. INTRODUCTION 

Millimeter-wave and infrared observations have long 
established that all stars form in molecular clouds; most 
of these form in Giant Molecular Clouds (GMCs). This 
conclusion has good theoretical support because Jeans 
instabilities with masses in the range of 0.1 - 100 M Q 
occur in regions with temperatures and densities found 
only in molecular clouds. Although good progress has 
been made in understanding how low-mass stars from, 
the physics of high-mass star and cluster formation re- 
mains elusive. Yet, nearly all of the information we have 
about present-day star formation in other galaxies comes 
from high mass stars. As a result, it would seem that un- 
derstanding star formation in other galaxies would be a 
daunting task. 

Despite our ignorance, it may not be necessary to know 
the details of star formation to understand how stars 
form on galactic scales. For example, a number of studies 
suggest that in the disks of normal galaxies, the efficiency 
of star formation — the fraction of the molecular mass 
of a galaxy turned into massive stars within 10 8 yr - 
shows an rms variation of only a factor of 2 when aver- 
aged over kiloparsec scales (e.g. iMurgia et al.l2002l etc.). 
This constancy suggests that understanding star forma- 
tion on kiloparsec or larger scales within galaxies may re- 
duce largely to a question of how molecular gas and thus 
how GMCs form. That is, once GMCs form, the average 
rate of star formation is determined to within a factor of 
2 as long as the IMF is relatively constant. Thus, finding 
a prescription for how GMCs form in a galaxy, together 
with a robust value of the star formation efficiency, may 
be all that is needed to determine the global star forma- 
tion history of a galaxy. Work of this sort has been pio- 
neered by Kennicutt (1983, 1989, 1998), but his inclusion 
of atomic gas, which is inert to star formation, is unsatis- 
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factory, since no known stars form from atom ic gas. Re- 
cently, however, iKrumholz fc McKed l)2005(l . have pro- 
vid ed a treatment th at justifies including atomic gas in 

the IKennicuttl l)1998fl star for mation prescript ion 

Following work by IWong fc Blitzl (|2002fl . 
iBlitz fc Rosolowskvl l)2004 hereafter BR04) investi- 
gated whether the mean ratio of molecular to atomic 
hydrogen, i? mo /, at a given radius in a spiral galaxy 
is determined primarily by a single parameter: the 
interstellar gas pressure, P ex t- They showed that this 
hypothesis leads to a prediction that the transition 
radius, Rt, where R m oi = 1, should occur at a constant 
value of £*, the stellar surface density. For a sample of 
30 galaxies, the stellar surface density at the transition 
radius (£*,*) is constant to within 40%, consistent with 
the hypothesis that R mo i is a function (/) of P ex t alone, 
i.e . flni.ni = f(Pe .T. t)- In this paper we extend the work 
of IWong fc Blitz! l)2002fl with a pixel-by-pixel analysis 
using additional galaxies to find the functional form 
of f(Pext) over a much larger range of pressure. We 
show that rms deviations from the derived relation are 
no more than about a factor of 2 for a range of galaxy 
types and metallicities. Finally, we use the pressure- 
molecular fraction rel ation to d e rive a star-formation 
law that modifies the IKennicuttl l)1998fl star formation 
prescription, particularly at low molecular gas fractions. 

2. BACKGROUND 

BR04 estimate P ex t for disk galaxies from the mid- 
plane pressure in an infinite, two-fluid disk with locally 
isothermal stellar and gas layers. They assume that the 
gas scale height is much less than the stellar scale height, 
which is typical of disk galaxies; so that, to first order: 
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Here, E 9 is the total surface density of the gas, v g is the 
vertical velocity dispersion of the gas, is the midplane 
volume density of stars, and p g is the midplane volume 
density of gas. 
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In most galaxy disks, p* is much larger than p g when 
azimuthally averaged. For a self-gravitating stellar disk, 
the stellar surface density, E* = 2p*ft*, where ft* is the 
stellar scale height and ft* = (v* 2 /27rGp*)°' 5 - Thus, ne- 
glecting p g , Equation 1 becomes: 



P ext = 0.84(GE*) a5 E 



(2) 



Direct solution of the fluid equations by numerical inte- 
gration shows that this approximation is accurate to 10% 
for £* > 20 Mq pc~ 2 (where p* < p g ), which covers the 
range of stellar surface densities in this study. 

We choose to express the midplane pressure in the form 
of Equation [3 since there is good evidence that both ft* 
and v g are constant with radius in disk galaxies (See 
BR04 for references). Furthermore, because of the weak 
dependence of P ex t on ft*, and t he small variation o f ft* 
measured from galaxy-to-galaxy ijKreeel et al.ll2002^ . we 
expect variations of ft* to have little effect on P ex t- 

By assumption, 
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E(H 2 )/S(H I) = f(P ext ) 



(3) 



Thus, since (v g /y/K^) is approximately constant within 
galaxies, then 



E(H 2 )/E(H I) = / [P ex t (Eg, £*)] 



(4) 



We wish to find / from Eq. Elfor as many galaxies and 
galaxy types as possible. 

3. ANALYSIS 
3.1. Data 

We have collected a sample of 14 galaxies for which 
the atomic, molecular and stellar surface densities can 
be determined. The galaxies and their adopted orienta- 
tion parameters are listed in Table ^ We also list the 
original publication of the H I and CO maps. Most of 
the H I data are from the VLA and the most of the 
CO data are from the BIMA Survey of Nearby Galax- 
ies l|Helfer et al.M2003D . For all of the galaxies in Table 
Q] except IC10 and the Milky Way, we derive their stel- 
lar surface densities from the 2 um (if s -band) ima ge in 
the 2MASS Large Galaxy Atlas ijjarrett I et all2 003'). We 
augment this samp le with the 2MASS K s image used in 
iLerov et alJ l)2006T ) for IC10 and data from the literature 
for the Milky Way. 

3.2. Derivation of Physical Quantities 

We measure the variation of molecular gas fraction 
(Rmoi = Eh2/Ehi) as a function of midplane hydrostatic 
pressure on a pixel-by-pixel basis from 2-dimcnsional 
maps of these quantities for the 14 galaxies in Table 
We estimate the midplane pressure in each galactic disk 
using an alternative form of the expression of BR04: 
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We determine the surface densities of the atomic gas, 
the molecular gas, and the stars from the maps of 21- 
cm, CO(J = 1 — > 0), and i4T s -band emission respectively. 



First, we remove foreground stars from the K s image 
by isolating point sources and replacing them with the 
median sky brightness surrounding each source. We then 
convolve the maps to the coarsest resolution among the 
three images (usually the H I). The resulting resolution 
measured in the plane of the sky is given in Table Q 
We then flag all pixels in the resulting maps that have 
less than a 5a rms significance in any of the integrated- 
intensity maps and do not consider these points in our 
analysis. Finally, we resampled the images onto the same 
coordinate grid so that the pixels in different images can 
be compared directly. 

Using the resulting maps, we estimate the surface den- 
sities of the various galactic components. For the atomic 
gas, the surface density is: 
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where i is the inclination of the galactic disk. We have 
assumed a mean particle mass of 1.36 m# to account for 
the presence of helium and metals. Similarly, the surface 
density of molecular gas is 
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adopting a CO-to-H2 conversion factor of N(H.2)/ Ico = 
2 x 10 20 cm" 2 (K km s -1 )- 1 and a mean particle mass of 
1.36 niH per hydrogen nucleus. Finally, we determine the 
stellar surface density from integrated K s -b&nd. light. 
We adopt a stellar mass-to-light ratio of Mk/Lk — 
0.5Mq/ Lq which approximates a l arge variety of star for - 
mation histories to 20% accuracy ijBell fc de Jongl200H) . 
Then 



logi 



E* 



Mq PC" 



= -0.4^a- s + 9.62 + log 10 cos i, (8) 



where p,Ks is the surface brightness of the galaxy in mag- 
nitudes per square arcsecond. 

Calculating P ext requires knowledge of the gas ve- 
locity dispersion and the stellar scale height. For all 
galaxies, we assume a constant gas velocity dispersion 
of v g —8 km s _1 , whi ch is observed across the disks 
of several galaxies Ce.g. iShostak fc van der KruiUlTosl 
Pic kev et alJlTflfl?! IBurtonlll 97lHMa,lhotrall1 9 95) . We es- 
timate the stellar scale height ft* based on the relation- 
ship between the radial scale length of the stellar disks 
(R*) and the corres ponding stellar scale heights found by 
iKregel et al.l ([2002) . Fitting their data gives 

log 10 ft* = (-0.23 ± 0.05) + (0.8 ± 0.1) log 10 R*. (9) 

where ft* and i?* are measured in parsecs. To measure 
i?*, we azimuthally average the stellar surface density 
in annuli of constant galactocentric radius assuming the 
orientation parameters in Tabled We fit an exponential 
function to the resulting profile ignoring regions indicat- 
ing the presence of a bulge. Using the derived value of 
i?*, we then calculate the scale height using equation 
The derived stellar scale heights are reported in Table 
HI The estimate of the scale height is calculated prior 
to convolution to a common resolution. We assume that 
that the ionized gas contributes negligibl y to E „; in any 
event, its scale height is large dRevnoldsl lT9891. and its 
contribution to the midplane pressure should be small. 
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TABLE 1 
Galactic Parameters 



Galaxy Name Distance Inc. P.A. H I data CO data Resolution Stellar Scale 

Height 

(Mpc) (°) (°) reference reference (kpc) (kpc) 



NGC 598 


0.85 


56 


203 


1 


2 


0.26 


0.21 


NGC 3521 


7.2 


58 


164 


3 


4 


1.3 


0.23 


NGC 3627 


11.1 


63 


176 


5 


4 


1.3 


0.39 


NGC 4321 


16.1 


32 


154 


6 


4 


1.6 


0.50 


NGC 4414 


19.1 


55 


159 


7 


4 


1.6 


0.31 


NGC 4501 


16.0 


63 


140 


6 


8 


1.8 


0.41 


NGC 4736 


4.3 


35 


100 


9 


4 


0.31 


0.11 


NGC 5033 


18.6 


62 


170 


10 


4 


2.1 


0.25 


NGC 5055 


7.2 


55 


105 


11 


4 


0.45 


0.24 


NGC 5194 


7.7 


15 





12 


4 


0.30 


0.33 


NGC 5457 


7.1 


27 


10 


11 


4 


0.20 


0.39 


NGC 7331 


15.1 


62 


172 


13 


4 


0.45 


0.28 


IC 10 


0.95 


48 


132 


14 


15 


0.28 


0.30 



References. - (1) IDeul fc van der Hulstl il98?h; £2 1 IHever et all <2004h: (3 1 M. 
Thornley, pri vate communic ation; (41 iH elfer et alJ 12003): (5) IZhang et al.l 119931): (61 
IKnapen eta!] <1993T); (7) IThornlev fc Mundvl 11997b!) : (8 1 IWong et alJ 120041: (91 IBraud 
119951): (10) IThean etaTI ll9971): (11) IThornlev fe Mundvl 11997al): (121 IRots et al.l 119901) : 
(13) IRegan et alJ 120041) 7(14) IWilcots fe Miller! 119981) (15) ILerov et all 12006!) 
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Fig. 1. — Molecular gas fraction as a function of midplane hy- 
drostatic pressure for NGC 5194 (M51). The mean relationship is 
plotted in a thick gray line. The completeness limit arising from 
the molecular gas observations is plotted as a solid black line. The 
intersection of this line with the bottom of the observed distri- 
bution determines the truncation limit. The relationship between 
P ex t and R mo i is only calculated above this limit. 



3.3. Completeness 

We calculate R mo i as a function of pressure on a pixel- 
by-pixel basis for all pixels that have measured surface 
densities of H I, CO, and 2 /xm emission. Figure ^ is a 
plot of the correlation between these two quantities for 
NGC 5194 (M51). The two quantities are clearly corre- 
lated, but the different signal-to-noise ratios of the three 
data sets affect the distribution of points in the Figure. 
Of the three tracers, the CO emission has the lowest 
signal-to-noise, and is found at the fewest positions. It 
is thus the limiting factor in our analysis. 

To illustrate the effects of finite sensitivity, we plot the 
locus of points where Sh2 = 8 M pc -2 (the 5a sensitiv- 



ity limit), = 300 M Q pc~ 2 , and a range of Shi as the 
solid curve in Figure^ We choose = 300 M® pc~ 2 , so 
that the curve roughly coincides with the left boundary of 
the data. The curve is a good match to the shape of the 
left edge of the data in Figure At a given pressure, 
points that lie below the line have CO surface bright- 
nesses that are too weak to be detected; the solid line 
therefore represents the completeness limit for the obser- 
vations. Because we wish to bin the data in values of 
Pext, we have complete CO data only when the lowest 
measured values of R m oi for the entire sample lie above 
the completeness line. Thus all of the data to the right of 
the dotted line, the intersection of the completeness line 
with the lowest values of R m oi have complete CO data 
(i.e. complete measured values of Sna)- Plots of P ext 
vs. Rmoi for four other representative galaxies are shown 
in Figure [3 showing that the scatter of the individual 
points within a galaxy exhibits considerable variation. 

3.4. The Relationship Between P ex t and R m oi 

Since we wish to measure the molecular gas fraction 
as a function of P ext /k, we average R mo i in bins of 0.1 
dex in P ext /k. To avoid biases due to the finite sen- 
sitivity, the analysis is restricted to P ex t/k larger than 
Pmin/k, the value where the completeness curve inter- 
sects the bottom of the data distribution. Practically, 
we define P m in/k as the bin containing the minimum 
value of Rmoi in the sample, and only consider bins of 
Pext/k with this value and larger. Since the upper limits 
on the CO surface brightness affect both P ext and R m oi , 
points without identifiable CO emission provide little ad- 
ditional information, and we discard these data from our 
analysis. The uncertainty in logi? mo ; (5 log Rmoi) is de- 
termined by the standard deviation of points in each bin 
(flog Rmoi ) divided by the square root of the number of 
independent samples in each bin: 

Slog Rmoi = ° X °J* m ° l , (10) 
V N l n beam 

where N is the number of pixels contributing to each bin 
and n oeam is the number of pixels per independent beam. 
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Fig. 2. — Plots of Pext vs. R mo i for four galaxies showing a range of scatter on a pixel-by-pixel basis. For galaxies such as NGC 3521, 
the scatter is so small (15%) that finding the completeness limit is trivial. For galaxies such as NGC 7331, we must use the technique 
described in the text for estimating the completeness limits in the relationship. 



We plot the relationship between hydrostatic pressure 
and molecular gas fraction in Figure [3J The figure shows 
that a power law is a good representation of the func- 
tional form between the two quantities. The plots for 
the individual galaxies have a similar slope but seem to 
lie in two groups. We fit a linear relationship between 
\og(P ex t/k) and log i? mo z to estimate the parameters for 
a power-law relationship between the two variables: 



Rmol — 



Pext 



(ii) 



The derived values of the index a are reported in Table 
with uncertainties generated by bootstrapping. Pq is the 
external pressure in the ISM when the molecular fraction 
is unity. The results for all of the galaxies is given by: 



Ptnol — 



Pext/k 



(3.5 ±0.6) X 10 4 



0.92±0.07 



(12) 



The data for NGC 3627, NGC 4321, and NGC 4501 
lie significantly above the data for the remainder of the 
sample. Splitting the data into two groups, we find 
(Po/fc) = 7700 cm -3 K for these three galaxies but 



(Po/k) = 43000 cm~ 3 K for the remainder of the sample. 
However, the mean index, (a), for each of the subsets is 
indistinguishable from that of the sample as a whole. 
The difference between the two populations can be at- 
tributed to differences in the H I content of the disks. In 
Table [3 we also tabulate the median values of Ehi cos i 
used in the analysis. We note that (Shi cos i) is signifi- 
cantly lower for these three galaxies than for the remain- 
der of the systems. The low values of Shi cos i likely 
arise from ram pressure or tidal stripping in these galax- 
ies. The other galaxies in our sample appear to be field 
galaxies without significant tidal or ram pressure inter- 
actions 3 . NGC 3627 is in the Leo Triplet of galaxies and 
has had a recent interaction with NGC 3628, evidenced 
by the asymmetric H I distributio n in NGC 3627 an d 
the long tail of H I from NGC 3628 ijrlavnes et al.lll979j) . 
NGC 4321 and NGC 4501 are both members of the Virgo 
Clust er and are known to be H I poor from previous 
work ijGiovanelli fc HavnesHl983|) . The H I morphology 

3 NGC 5194 (M51a) is interacting with NGC 5195 (M51b), but 
the interaction shows no evidence of stripping neutral gas from 
NGC 5194. 
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Fig. 3. — Molecular gas fraction as a function of midplane hydrostatic pressure for 14 galaxies. All galaxies follow a similar scaling 
between P ex t and R mo i- The mean index of the relationship is 0.92 ± 0.07 (dashed line). There are two categories of galaxies, offset from 
each other by a factor of 5 in Pq . 



TABLE 2 

Results of Pressure Analysis 



Galaxy 


a 


P /k 


Scatter 3 - 


(Ehi cos i) 


Morphological 


Name 




(10 4 cm" 3 K) 




(M© P c- 2 ) 


class b 


MW 


1.64 ±0.11 


2.0 


0.09 


8.4 


Sb 


NGC 598 


0.87 ±0.14 


5.1 


0.03 


9.3 


SA(s)cd 


NGC 3521 


1.02 ±0.03 


7.1 


0.02 


16.8 


SAB(rs)bc 


NGC 3627 


0.81 ±0.03 


0.1 


0.10 


4.3 


SAB(s)b 


NGC 4321 


0.84 ± 0.03 


0.7 


0.06 


6.7 


SAB(s)bc 


NGC 4414 


0.89 ±0.02 


4.6 


0.02 


14.4 


SA(rs)c? 


NGC 4501 


1.07 ±0.08 


1.2 


0.13 


4.2 


SA(rs)b 


NGC 4736 


0.93 ±0.04 


6.5 


0.09 


11.1 


(R)SA(r)ab 


NGC 5033 


0.76 ± 0.03 


3.0 


0.05 


11.5 


SA(s)c 


NGC 5055 


0.73 ± 0.02 


2.S 


0.03 


11.5 


SA(rs)bc 


NGC 5194 


1.00 ±0.05 


3.0 


0.07 


13.2 


SA(s)bc pec 


NGC 5457 


0.58 ± 0.32 


2.1 


0.09 


16.6 


SAB(rs)cd 


NGC 7331 


1.01 ±0.06 


5.1 


0.05 


14.5 


SA(s)b 


IC 10 


0.73 ±0.14 


5.6 


0.10 


6.4 


dlrr IV/BCD 


Mean 


0.92 ±0.07 


3.5 


0.06 


10.6 


Sbc 


Mean (Non-interacting) c 


0.92 ±0.10 


4.3 


0.05 


12.2 


Sbc 


Combined Data 


0.94 ±0.02 


4.5 


0.14 


9.9 





a Dcfined as the standard deviation of the residuals: A log R mo i = log R mo i — a(log P e xt — log Po) 

b From RC3 Ide Vaucouleurs et alJITTmilh 

c Excludes data from NGC 3627, NGC 4321 and NGC 4501. 
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Fig. 4. — Molecular gas fraction as a function of midplane 
hydrostatic pressure for all data in the sample. Each point repre- 
sents an independent measurement of Sjji, Sh2 an d S*. A linear 
fit to log_R mo ; vs. log Pext/k is shown as a thick gray line. The 
data for the interacting galaxies (red, NGC 3627, NGC 4321, NGC 
4501) are excluded from the fit. The results are indistinguishable 
from determining the relationship between R mo i and P ex t for each 
galaxy individually. 



of NGC 4501 shows ample evidence that the galaxy is 
undergoing ram pressure stripping in the high-pressure 
cluster environment. NGC 4321 is less obviously dis- 
turbed, but still H I poor. The three galaxies do not 
appear to be discrepant in any other way. 

If we remove these galaxies from the mean relationship, 
we obtain: 



R>i 



Pext/ k 



(4.3 ±0.6) x 10 4 



0.92±0.10 



(13) 



which changes Pq by only 20%. This may be a more rep- 
resentative relationship than Equation 1121 because the 
discrepant galaxies may be subject to pressure in addi- 
tion to the hydrostatic pressure, which would move them 
to the right in Figure |21 We also fit a power law to the 
combined data for the non-interacting galaxies, rather 
than fitting each galaxy and then averaging (Figure 
This weights the fits by the total number of data at a 
given P ex t- Because the nearby, H I-rich galaxies (IC10, 
M33) have many more independent samples, this fit em- 
phasizes galaxies with low values of R m oi- Nonetheless, 
the resulting fit (Combined Data in Table |2J) is indistin- 
guishable from the average of the individual fits to the 
non-interacting data. 

It is worth noting that although we derive a relation- 
ship involving the gas pressure, because we take the ve- 
locity dispersion of the gas as constant, and because the 
turbulent pressure, pv 2 , provides most of the support 
of the gas layer, our expression is essentially equivalent 
to Rmoi being a function of the midplane gas density. 
We cannot, with our data alone, decide between these 
two possibilities. But using the equivalent density pa- 
rameterization yields a true Schmidt law, wher e the star 
forma tion rate is dependent on ISM density ijSchmidtJ 
1959). 



3.5. Comparison with Previous Work 
iWong fe Blit3 p00l found, for a sample of seven 
galaxies, that R mo i oc P ex t°' 8 over one and a half or- 
ders of magnitude in pressure. For a given R mo i, Pext 
ranges over an order of magnitude suggesting that the 
dispersion in the pressure is about a fact or of 2 - 3. We 
have m ade several improvements to the IWong fc Blitzl 
( 2002) analysis as follows. Rather than taking azimuthal 
averages, we h ave calculated R mn i a nd P ex t on a pixel- 
by-pixel basis. IWong fc Blitzl i|2002fl used an expression 
for the pressure by ? that requires knowledge of the stel- 
lar velocity dispersion in disks, a quantity that is difficult 
to measure, and one that is expected to vary significantly 
as a function of radius. In this paper, we use a different, 
but equivalent expression that relates the pressure to the 
stellar surface density and stellar scale height (BR04). 
The surface density is observable; the scale height varies 
little within galaxies and it s value can be infe rred from 
the observable scale length l)Kregel et al.l l2002. Equation 
B y including addition al galaxies, we are able to ex- 
tend IWong fc Blit3 <|2002fl to a total of three orders of 
magnitude and to include lower mass and lower metal- 
licity galaxies. Our method also significantly reduces the 
scatter among the galaxies, measured by the dispersion 
of the values of Po/k in Table |3 

We note that over the range of galactic radius for M33 
included in our samp le, there is a decrease in the metal- 
licity by a factor of 6 (H enrv fc Howard 19951) . Although 
the range of pressure sampled in M33 is small, there is no 
systematic offset from the mean relation shown in Figure 
13 with metallicity. The same is true f or IC10, which has 
a metallicity 1/5 solar llGarnettlll990l). 

On theoretical grounds. lElmegreenl l)1993f) derives a 
much steeper dependence of R mo i on pressure than we 
do, but he also includes a dependence o n the emissiv- 
ity of the dissociating radiation field, j. IWong fe Blitzl 
(2002) found that if one assumes that j oc T, gas oc E*, 
Elmegreen's (1993) expression is equivalent to f mo i 



Ef^/Eg P e . 



1.7 



However, IWong fe Blit3 (|2002|) did 



not take into account that c*, the velocity dispersion of 
the stars, has a radial dependence on E*, which is ap- 
proximately oc E* ' 5 . By ignoring the gravity of the gas 
as we have done in equation [21 and including the depen- 
dence of c* on E+, Elmegreen's expression for P, has the 
same functional dependencies as ours. In the regime of 
low fmoi ! where Elmegreen's treatment is meant to be 
appropriate, f mo i ~ R m oi- Although j probably does not 
have a simple power law dependence on S ffas , if it can 
be approximated as j oc E g , then Elmegreen's (1993) ex- 
pression for f mo i has a power law dependence quite close 
to what we find observationally. 

We note further that although one expects R mo i to 
have a dependence on j on physical grounds, the small 
scatter of the lines in Figure [3] implies that empirically, 
Rmoi behaves as if it depends on pressure only. The 
discussion above suggests that the dependence on j is 
subsumed within the empirical relationship. 

3.6. Notes on Individual Galaxies 

The Milky Way - For the Galaxy, we use the values of 
iDamel l|1993fi for E H i and E H 2 at R ga i > 3 kpc. We as- 
sume that the local stellar surface density in the vicinity 
of the Sun is 35 M Q pc -2 and that the stars are dis- 
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tributed in an exponenti al disk with scale length 3 kpc 
Ipehnen fc Binnevlll998h with R = 8.0 kpc. We scale 
Eh2 by 1.36 to account for the presence of helium and 
by (2/2.3) to account for the different conversion factor 
adopted in that study. Finally, we scale Ehi by 1.5 to ac- 
count for the presence of high latitu de gas, based on th e 
results beyond the Solar circle from lLevine et all {2006). 

NGC 598 (M33) - NGC 598 is the closest extragalac- 
tic source in our sample. The disk of NGC 598 is stable 
accor ding to the Toomre criterion ijMartin fc Kennicutti 
|2001|) . making it a notable test case for for comparing 
the role of gravitational instability with that of pressure 
in regulating the ISM and star formation. We gener- 
ated a molecular gas map from a masked integrated- 
intensit y map from the FCRAO 14-m survey of the 
galaxy l)Hever et all2004l) . Apart from the Milky Way, it 
is the only other disk galaxy in the sample with sufficient 
surface brightness sensitivity in the CO map to measure 
Rmoi where R m oi < 1- The de rived index of a = 0.87 
can be compared to the work of iHever et al.l (|2004fl who 
show that R mo i oc P^ t i~ x where j is the mean radia- 
tion field. Using j cx E*, we reproduce their index to the 
pressure relationship. Performing a similar analysis for 
the remainder of our sample finds indices ranging from 
1.4 to 2.3. 

NGC 3627 - This galaxy is a molecule rich member 
of the Leo Triplet and the ISM of the galaxy has been 
affected by interaction with the remainder of the group. 
NGC 3628, a nother member of t he triplet, is also very 
molecule rich ijYoung et al.lll99"5l) . The interactions be- 
tween the three galaxies have stripped a long plume of 
material from NGC 3628 and have also likely depleted 
the H I in the disk of NGC 3627. 

NGC 4321 (M100) - This galaxy is a member of the 
Virgo cluster, but is not obviously undergoing ram pres- 
sure stripping at the present time. Azimuthally averag- 
ing the H I data shows a slight downturn at large radii, 
and the median surface density for the H I is relatively 
high (6.7 Mq pc -2 ) compared to the other interacting 
galaxies in our sample (NGC 3627 and NGC 4501). Of 
the three galaxies undergoing interactions, this galaxy 
appears to be the least affected by interaction. How- 
ever, it is noticeably deficient in H I with a total atomic 
gas mass a factor of three lower than comparable field 
galaxies ijKennev fc Youndll989|) . 

NGC 4501 (M88) - We used archival BIMA data from 
NGC 4501 to generate the CO map of NGC 4501 used 
in this study. The data were originally presented in 
iWong fc Blits3 <|2002). NGC 4501 is a member of the 
Virgo Cluster and its H I morphology appears to be no- 
ticeably affected by the cluster environment. It is sepa- 
rated by ~ 500 kpc on the plane of the sky from M87; and 
the H I surface density is asymmetrical, a pparently piled 
up on the side of the galaxy closest to M87 ijCavatte et al.l 
1990). The cluster environment appears to be responsi- 
ble for the depletion of the atomic gas which is a factor of 
3 lower than what wou ld be expected from a fi eld spiral 
of similar stellar mass l|Kennev fc Younell 1989ft . 

NGC 4736 (M94) - This galaxy has the earliest mor- 
phological type in our sample and sizable amounts of 
neutral gas are found within the bulge of the galaxy. 
The pressure formulae used in this paper is appropriate 
for the disks of galaxies, but should be an upper limit 
for the bulges of galaxies, since the thickness of bulges 



is larger than disks. Applying the pressure formulae for 
disks to NGC 4736 may account for the high value of Pq 
in the system. 

NGC 5194 (M51a) - NGC 5194 is a classic example 
of grand-design spiral structure in galaxies. The high 
resolution observations of NGC 5194 facilitate the mea- 
surement of P ext and R mo i both in and out of the spiral 
arms seen in the ISM. In both the arms and the interarm 
regions, P ext appears to be a good predictor of Rmoi- 

NGC 5457 (M101) - The overlap between H I and 
CO detections spans only a very small portion of the 
galaxy. The H I surface density drops markedly in the 
center of the galaxy where the molecular gas is concen- 
trated ijWong fc Blitzll2002|) . In the small overlap region, 
however, the data from NGC 5457 agree well with the re- 
mainder of sample. 

IC 10 - We include IC 10 in our sample to investi- 
gate the Rmoi vs. P ex t relationship in a low metallicity 
environment. The gas geometry of the system is quite 
complex and not apparently confined to the stellar disk, 
whic h may affect the ap plicability of Equation |3J How- 
ever, [Lero^^OIQ <|2006|) noted that the molecular gas in 
the system tended to be found on filaments of atomic gas 
that were coincident with the stellar disk of the galaxy. 
We include the pressure data from their work in this 
study. They assume a characteristic scale height of the 
stellar component of 300 pc. We correct the surface den- 
sities of gas by the inclination of the stellar disk assuming 
that the gas spatially coincident with the stellar disk ac- 
tually lies within the disk. 

3.7. Systematic Effects 

In this section, we discuss how variations in our as- 
sumptions and interpretations affect the results. 

Since Eh2 appears in the expressions for both R mo i and 
Pexti it is possible that some of the observed correlation 
in Figure |21 arises simply because we are correlating CO 
emission with itself. As we can see from Figure 01 how- 
ever, most of our individual data points come from the 
low pressure regime (P < Pq), where no such correlation 
is expected. Thus, this concern is most important where 
Rmoi 3> 1- In that case, the Equation |21 as plotted in 
Figure El reduces to: 

(X 

logS H 2 - logSffl = a logS H 2 + g lo S s * + const. (14) 

where a is the quantity we are trying to determine (as- 
suming that Rmoi is proportional to P a ). Since Shi 
shows little variation compared to either Eh2 or E*, and 
because log Eh2 dominates the left hand side of Equation 
ITU for large R m oi , we can absorb log Ehi into the con- 
stant. Because log Eh2 appears on both sides of Equation 
1141 at least some of the correlation at R mo i ^ 5 is due to 
correlating log E H 2 with itself. The slope of the P — R mo i 
relationship is given by: 

a _ log E H 2 + const. 

a ~ logE H 2 + (logE*)/2 1 ' 

As long as the range of values of (logE*)/2 are compa- 
rable to the range of values of log Eh2 > as is the case for 
most of the data presented in Figure|21 both terms in the 
denominator will determine the slope, and thus the value 
of a. Thus, even for large values of R m oi, where some 
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of the correlation in the P — R mo i relation is produced 
by Er2 appearing on both axes, the data still provide 
information about the functional form of Equation |3J 

To estimate the amount of correlation that arises 
purely from the mathematical construction of our anal- 
ysis, we recalculate R mo i and P ext after scrambling the 
original measurements of £hi,£h~2 ano - This ran- 
domization retains the original data but the surface den- 
sity measurements for each disk component no longer 
correspond to the values of the other components. Us- 
ing the scrambled data for R mo i and P ext we do indeed 
find a significant correlation in the randomized data, par- 
ticularly evident at R m oi ^ 2. The power-law index 
for the scrambled data is a — 1.5 ± 0.02, significantly 
greater than the derived relationship for the real data. 
We conclude that we are observing a physically signifi- 
cant relationship although we recognize that the data for 
Rmoi ^ 2 would not be compelling on their own. When 
used in conjunction with the data for low values of R mo i, 
however, the molecule-rich regions of galaxies extend the 
trend into high-pressure regions. 

Another concern is that variations in the CO-to-fi2 
conversion factor may substantially alter our results. 
However, the only regions included in this sample that 
show good evidence for variations in the CO-to-H2 con- 
version factor that are also included in this sample are the 
mole cule rich regions in th e central ~ 200 pc of galaxies 
(e.g. iMeier fc Turnerl l2004). In these regions, R m oi 3* 2 
and the observed correlation is largely the product of cor- 
relating CO emission with itself. As a result, changes in 
the CO-to-H2 conversion factor will alter R mo i and P ex t 
by roughly the same amount and the data will likely 
remain consistent with the observed trend. Thus, varia- 
tions in the conversion factor are unlikely to substantially 
affect our results. 

The projected resolution at the distance of the galaxies 
varies by nearly an order of magnitude across the sam- 
ple. There is no significant correlation between physical 
resolution and the index of the derived relationship or 
the degree of scatter around the mean trend. To assess 
the effect of the physical resolution on our results, we 
repeat the analysis, convolving the data to a constant 
physical resolution of 1.8 kpc. This degrades the resolu- 
tion to a constant scale for all galaxies except for NGC 
5033. We find no significant change in the results as a 
result of smoothing to this resolution. We conclude that 
resolution effects do not affect the results significantly. 

The assumption of a constant velocity d isper- 
sion is well-justified ifvan der Kruit fc Shostakl IT982I: 
[ Shostak fc van der Kruirill984Hvan der Kruit fc ShostaU 
l!984Hvan Zee fc Br vantl 1999ft . but there are ~ 30% vari- 
ations in velocity dispersion within and among galaxies. 
These differences may account for the offsets and scatter 
seen in Figure 

4. DISCUSSION 

In nearby galaxies where the dominant state of the gas 
is atomic, the atomic gas occurs pr i marily in filaments 
llKim et alJ Il998t lEngargiola et alJ l2003t iLerov et all 
2006), and the formation of the filaments appears to 
prece de the formation of the molecular gas l)Blitz et al.l 
•JOOfi;. From §3 and BR04, we have shown that hydro- 
static pressure alone is sufficient to determine what frac- 
tion of the neutral gas at any location in a galaxy is 



molecular, and we have determined the quantitative re- 
lationship between pressure and molecular gas fraction. 
We have not been able to show, however, how the molec- 
ular gas collects itself into GMCs, or even if the atomic 
gas first collects itself into GMC-sized objects and then 
becomes molecular. This detail may not be important 
for determining the star formation rates in galaxies on 
kiloparsec scales, since the rate at which stars form from 
the molecular gas ap pears to be well dete rmined on scales 
larger than a GMC IjWong fc Blitzll2002j) . but is relevant 
to issues related to the stability of the ISM. 

In all of the galaxies we have examined in this paper, 
pressure is measured on a scale large compared to the size 
of a GMC, because of either the relatively poor resolu- 
tion of the H I observations, or because the galaxies are 
relatively distant. We expect that Equation El would 
break down on a scale comparable to that of a typical 
GMC because the pressure within a GMC is enhanced by 
its self-gravity and is no longer in pressure equilibrium 
with its environment. We would guess that the scale on 
which the relation is valid is about equal to twice a pres- 
sure scale height (~ 400 pc in the solar vicinity of the 
Milky Way) . The pressure will be in equilibrium on such 
a vertical scale, on average, and is unlikely to vary by 
very much on a similar scale in the galactic plane. 

4.1. Applicability to Other Kinds of Galaxies 

Our sample of galaxies is limited to normal spirals and 
some dwarf galaxies. It contains weak Seyferts and galax- 
ies with low level starbursts and galaxies with metallici- 
ties that range from solar values to values as low as 1/6 
solar (IC 10). It does not include elliptical galaxies, and 
SOs, a number of which have been sh own to be molecule 
rich (|Youndl2002llYoung et aTlll995ft . nor does it include 
very low metallicity dwarfs such as I Zw 18, luminous 
starbursts, strong Seyferts or merging systems. Equa- 
tion El may not be applicable in these systems, but the 
pressure in the ISM can be estimated from other means. 
Equation 1131 uses an equilibrium approximation of disk 
structure to estimate pv g 2 , a quantity that can be deter- 
mined equally well using a variety of methods. In ellipti- 
cal galaxies, the pressure may be measured directly from 
density and temperature information inferred from x-ray 
data. In the bulges of disk galaxies, other constraints 
from stellar velocity dispersions and scale heights can be 
used to estimate a more accurate stellar mass distribu- 
tion. Programs are underway to to probe the relationship 
between molecular gas fraction and ISM pressure in an 
array of such systems. The Large and Small Magellanic 
Clouds will also be good testbeds for determining the 
general applicability of the pressure-molecular fraction 
relation. In these two systems with irregular morpholo- 
gies, we can measure the relationship between molecular 
gas fraction and pressure at high resolution. 

Galaxies with low metallicities are interesting regions 
in which to study the molecular fraction of the ISM. In 
these systems, the flux of dissociating UV photons tends 
to be higher than in the disk galaxies that we studied in 
this sample. In conjunction with the volume density of 
hydrogen, the UV flux establishes the equilibrium molec- 
ular fraction in the ISM of the neutral gas. Hence, we 
suspect systems with low metallicity and high UV flux 
(such as the SMC) may be significantly offset from the 
Rmoi ~ Pext relationship seen for disk galaxies with ap- 
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proximately solar metallicity. 

4.2. Molecular Gas and Star Formation 

IKennicuttl l|1998l hereafter K98) has argued that the 
surface density of the rate of star formation Ssfr is pro- 
portional to the sum of the atomic and molecular gas 
surface density, S g over a wide range of galaxy luminos- 
ity. This relationship is puzzling because stars form only 
from molecular gas. K98 cites the relative tightness of 
his plot of Ssfr vs. S ff compared to his plot of Ssfr vs. 
Sjj2- However, when averaged over kiloparsec scales, Shi 
tends to saturate at a value of ~10 M m p_ c ~ 2 ^ ° r about 
1.3 x lO 21 cm" 2 (e.g. iCavatte et al1IT9MlWong fc Blitzl 
2002), the relationship between Ssfr and Shi in K98 is 
less a correlation than a plot showing the saturation of 
t he H I line. 

iWong fc BliU l|2002l hereafter WB02), on the other 
hand, have shown that for 7 well-studied gala xies in the 
BIMA SONG CO survey, l|Helfer et a,U l200.U the pro- 
portionality found by K98 does indeed hold, but for the 
molecular gas only. That is, K98 finds that Ssfr = 
0.16 S 4 - 4 , where Ssfr is measured in M pc~ 2 Gyr -1 
and S g is measured in M Q pc~ 2 . WB02 find that 
Ssfr = 0.13 S H 4 . These relations use values of Xco 
that differ in whether or not they include a correction 
for helium. If helium is included in the WB02 relation, 
then their coefficient is 0.18 rather than 0.13, which dif- 
fers from the K98 value by only 13%. This is not entirely 
surprising, since WB02 use K98's calibration of Ha line 
st rength vs. sta r form ation rate. 

Mureia et al.l l)2002|) developed a star formation for- 
mulation based on radio continuum measurements cal- 
ibrated to the supernova rate in a number of galaxies 
ijCondonl 11992). This relation, like WB02, relates the 
star formation rate to Sjj2 and finds a power-law index 
of 1.3, indistinguishable from 1.4 within the uncertain- 
ties. The IMurgia et al.l l|2002|) coefficient of proportion- 
ality is a bit higher: 0.26 rather than 0. 18, but u s es the 
Mille r fe Scalol ((1979^ IMF rath er than alSalpeterl (fl955h 
IMF between 0.1 and 100 M . IMurgia et all <|2002t) do 
not give a lower mass limit, but the differences between 
the two IMFs could make up much of the discrepancy in 
the coefficient. 

The K98 relation is based primarily on global aver- 
ages, and the relation of WB02 is based on pixel-by-pixel 
measurements convolved to the same linear resolution. 
As a practical matter, the differences between these for- 
mulations are unimportant for regions of galaxies where 
the ISM is dominated by molecular gas, such as galac- 
tic centers, as well as most of the regions of galaxies over 
which BIMA SONG has enough sensitivity to detect CO. 
However, throughout the entire disk of th e Milky Way , 
Shi > Sh2 (except for the central 300 pc: | Dam e 1993), 
as is also the case in M33 lIHever et al.l l2004) and most 
of M31 ijNieten et al .1 12006). In general, the interstellar 
medium in dwarf galaxies is dominated by H I rather 
than m olecular gas, and it is i n M33 and in dwarf galax- 
ies that iMartin fc Kennicuttl l)2001fl find that their star 
formation prescription breaks down. Observations of H I 
rich dwarf galaxies should be a good test of whether the 
K98 or WB02 relation is better at predicting the star 
formation rate. On e recent study has been completed by 
iHever et al~l (|2004j) on the Kennicutt-Schmidt law in M33 



where they find Ssf r. oc S^ 6 in good ag reement with the 
work of WB02 and IMurgia et al J (|2002T . When they in- 
clude atomic gas in their analysis, they find Ssfr oc S h '|, 
differing dramatically from the results of K98. 

4.2.1. CO, HCN, FIR and Star Formation 

In her ground breaking work on the re lationship be- 
tween dense gas and star formation, lLadal l)1992fl showed 
that in the Orion B molecular cloud (L1630), nearly all 
of the star formation is associated with the H2 traced by 
the CS molecule, a tracer of dense molecular gas. Most 
of the molecular gas in L1630 traced by CO, a moderate 
density tracer, is inert and does not take part in the star 
formation process. Why, then, is there a good relation- 
ship between Sco and star formation if the molecular gas 
traced by CO h as little direct be a ring o n star formation? 

The work of IGao fc Solomon! ll200l hereafter GS04) 
and of IWu et all II2005D suggests an answer. GS04 show 
that HCN is a linear tracer of the star formation rate in 
luminous and ultraluminous galaxies while CO is not. In 
their work, the star formation rate is determine d from 
the far IR luminosity. The connection to Lada's (1992) 
work is that both HCN and CS are good tracers of dense 
gas with n(H2) ~ 10 4 ~ 5 , but H CN is somewhat more 
luminous, on average, in galaxies ijHelfer fc Blitzlll997bL 
and references therein) . The GS04 plots of the correla- 
tion between HCN and CO with IR luminosity show what 
appear to be comparably good fits to the data, though 
the fit using CO is not linear, while the fit using HCN is. 
The correlation coefficients for both fits are quite high, 
but the fit using HCN is marginally better than for CO 
(R = 0.94 vs. 0.88), and the scatter somewhat less (0.25 
dex vs. 0.35 dex). In luminous and ultra-luminous in- 
frared galaxies, GS04 find that the HCN luminosity is 
related to the CO luminosity: Lhcn = 1.6 x 10~ 5 ^co 8 
f or L measured in K km s _1 pc 2 . 

IWu et akl i|200"> ■ show that the linearity between the 
HCN and IR luminosity continues down to the scale of 
individual dense cores in Galactic molecular clouds, im- 
plying that the HCN - IR relationship for starbursts con- 
tinues down to the level of star-forming clumps. The lin- 
earity holds over a total range of luminosity of 7 - 8 or- 
ders of magnitude. Correlations using luminosity usually 
look remarkably good because one is correlating distance 
(squared) with itself on both axes, making the correla- 
tion much more striking than it actually is. A better 
comparison is between distance independent quantities 
such as surface brightness. Nevertheless, there seems to 
be good evidence that the linearity of IR luminosity with 
HCN line strength is quite good for a range of objects 
from individual Galactic GMCs to ULIRGs. 

4.2.2. A Pressure-Based Prescription for Star 
Formation on Galactic Scales 

We propose an alternative, pressure-based prescrip- 
tion for galaxy-scale star formation centering on three 
empirically-derived relationships. First, we assume that 
the FIR luminosity is a good tracer of the massive star 
formation in a galaxy because the FIR emission is pro- 
portional to the number of Lyman continuum photons, 
nearly all of which are absorbed by dust, From GS04 
(and references therein), 

M SF R«2 x 10- 10 (L FIR /L Q ) M q yr" 1 , and (16) 
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M, 



dense 



= 14 



L 



HCN 



K km s 1 pc 2 



Mr. 



(17) 



Next, we assume that all star formation takes place in 
dense molecular gas which has a constant star formation 
efficiency. Expressing this in a form comparable to that 
of K98 using the above notation gives: 

E SFR = eE H2 (HCN) (18) 

GS04 and lWu et all ll27)05ft measure e « 10 - 13 Gyr^ 1 . 
Expressing Equation ^] in terms of the total molecular 
gas traced by CO gives: 



Esfr. = eEh 2 (CO) 



E H2 (HCN) 



E H2 (CO) 

which is related to the total gas content: 

'E H2 (HCN) 



J SFR 



gjmol 



E H2 (CO) 



(19) 



(20) 



where f mol = E H2 /E g = R mo i/{R mo i + 1) and R mo[ = 
Eh 2 /Ehi- Third, we assume that fmol is determined from 
the pressure and the total amount of gas present, and can 
be calculated from Equations llllandlT2*l. We can rewrite 
Equation m 



fmol 



{Rmol + 1) 



Po 



(21) 



For an integrated intensity, J, expressed in units of 
K km s _1 and E measured in M Q pc -2 , Eh 2 (HCN) — 
14 J(HCN) (GS04), and E H2 (CO) = 4.4 /(CO) for our 
adopted CO-to-H 2 conversion factor ( H3.2| l corrected for 
helium. Thus, 



S H2 (HCN) 
E H2 (CO) 



3.2 



/(HCN) 
/(CO) 



(22) 



Unlike theoretic al explanations of the K ennicutt- 
Schmidt law (e.g. iKrumholz fc McKed feOOS). this for- 
mulation separates the empirical star formation relation- 
ship of K98 into three constituent empirical relationships. 
The ability to reproduce the K98 results should not be 
surprising - all three of the constituent relationships are 
based on galaxies similar to those in K98. A consistent 
constant of proportionality is an indication of accurate 
calibration of the various quantities involved. What is 
notable about this formulation is the ability to charac- 
terize star formation on a galactic scale using the star 
formation efficiency of dense gas and the role of pressure 
in determining both the fraction of gas that is molecu- 
lar and the fraction of that molecular gas that is dense 
enough to form stars. 

The Low Pressure Regime 

Except for the central 300 pc of the Milky Way, and the 
inner regions of M31, the mean H I surface density ex- 
ceeds that of the H 2 when averaged over several hundred 
pc everywhere in each galaxy. Thus, most of the molec- 
ular gas in these galaxies occurs in the regime where 
Po > P ex t- For most GMCs in Local Group galaxies, 
the surface density of GMCs as traced by C O is roughly 
const ant with a value of ~ 100 M Q pc -2 ifBlitz et al.l 
2006), so the fraction of dense gas and consequently, 
the star formation rate of individu al GMCs, will also 
be roughly constant, as is observed ijMoonev fc Solomon! 



Il988t IScovTlle fc Goodl ll989T) . Variations in Esfr over 
kiloparsec scales thus result from changes in the frac- 
tion of neutral gas bound up in GMCs (oc fmol)- That 
is, averaging emission on large enough scales (> 200 pc) 
over which the surface filling fraction of molecular gas is 
small, pushes the Kennicutt-Schmidt law into the regime 
of small values of Er 2 since averaging GMCs with atomic 
gas that has no star formation will geometrically dilute 
both Esfr and Ejj 2 (CO) by a comparable, though not 
identical, factor. 

In th is regime, we may tak e Equation 1221 to be con- 
stant. iHelfer fc Blit3 l|1997aft derive a value of 0.026 
±0.008 for the ratio of integrated intensities averaged 
over the inner disk of the Milky Way. Similarly, GS04 
find L(HCN)/L(CO) = 0.039 for their low-luminosity 
galaxies. We adopt the average of these two values as a 
characteristic line ratio which implies that the ratio in 
Equation [53 is 0.1. Using this and taking P ext <C Po, 
Equation [201 becomes 



J SFR 



0.1 eE fi 



Po 



M Q pc -2 Gyr" 



(23) 



where P ex t is given by Equation [5] and Po/k = 4.3 x 10 4 
cm K. This expression will, in general, give smaller 
values of Esfr than K98. 

The High Pressure Regime 
For galaxies where R m oi ^> 1 over large scales in the 
ISM, the situation changes. In this case, P ex t > Po, 
and fmol ~ 1. In at least part of this regime, GMCs 
do not have a constant surface density, and one cannot 



take Equation [53 to be constant. For /, 



mol 



> 



0.8, E, 



Eh 2 (CO). Then the observed Schmidt-law relationship 
in th e molecular gas (Esfr cx Eh 2 1,4 : iWong fc Blitzj 
I2002t iMurgia et alll2002t iHever et all 12004ft reduces to 
the prescription of K98, Esfr = 0.16 S ff 14 . 

An alternate expression can be derived from the con- 
side ration tha t stars form only from dense gas. GS04 
and lWu et al. (2005j) propose a star formation prescrip- 
tion based on the dense gas content of a region: 



M* = (1.1 ± 0.2) x 10- 8 M denS e M Q yr" 



(24) 



where Mdense is the mass of dense gas including helium 
measured in solar masses. Using the scaling between 
^hcn and Lqq used in the GS04 study, this suggests 
the star formation rate can be deduced from the total 
molecular gas mass of the galaxy: 



M* = (0.77 ± 0.07) 



/ M] 



H2 



V.10 9 m q 



1.44 



M Q yr" 



(25) 



In the high-pressure regime, Mh 2 « M g so M+ oc Mg A . 
While the index of 1.4 is reminiscent of that found in K98, 
the two expressions are qualitatively different. Equation 
[2*B1 deals with the global properties of a system (mass) 
rather than the local properties (surface density) used in 
K98. 

To what degree are the prescriptions, Eauations l23l aiid 
1251 both locally and globally applicable? Both equa- 
tions use relations that are globally determined (Equa- 
tions El and ^J, so they will both be applicable globally. 
But both equations also use locally determined relations 
fEauations 1121 and I20f> . As long as care is taken to nor- 
malize the global relations over the same area, the basic 
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relations should be applicable locally as well (over ~ kpc 
scales). Equation[23is determined from CO observations 
that include luminous IR galaxies (LIRGS - GS04), so 
should be applicable to these systems as well. However, 
because a scaling has not been established for ULIRGS, 
care should be taken when applying our relations to these 
systems. We use these relations both locally and globally 
m Section |4~Q below. 

4.2.3. Comparison with Observations 

In this section, we compare the results of our star for- 
mation prescription to observed estimates of the star for- 
mation rate in various galaxies. 

The Milky Way - We estimate the star formation 
rate in the Milky Way using Equation |201 We ex- 
press fmoi in terms of the pressure (Equation 1211) ; be- 
cause the disk of the Galaxy is entirely in the low 
pressure regime, we assume a co nstant dense gas frac- 
tion: E H 2(HCN)/E H 2(CO) = 0.1 ijHelfer fc Blitdll997at 
GS04). This gives an expression for Esfr as a func- 
tion of Galactic radius. We plot the prediction in Fig- 
ure El The surface density of star formation expected 
from the K98 prescription is shown as a dashed curve, 

We also com- 



J SFR 



= 0.16 Eg 1 ' 4 M Q pc~ 2 Gyr~ 



pare the results to the observed data for H II region s 
ilGuesten fc Mezgerl 11982ft. nulsars lILvrie fit, al] [l985), 
and s upernova remnants llGuibert et alJll978j) . as plot- 
ted bv lPrantzos fc Aubertl l|l995f) . normalizing to a local 
star formation rate: Esfr(-Ro) = 1-9 Mq P c_2 Gyr -1 
ijKroupa Ill995j) . In general, the pressure-based predic- 
tion agrees reasonably well with the observed data. In 
contrast, the K98 prediction gives good agreement at 
the peak of the molecular ring (R ga i ~ 4 kpc). but 
does not predict the shape or the amplitude of the curve 
well at larger R ga i; it overpredicts the SFR, as expected. 
Neither does particularly well inside the peak of the 
molecular ring, but that may be due to uncertain kin- 
matics induced by the bar. In the Solar neighborhood, 
we predict Esfr = 2.2 M Q pc -2 Gyr -1 compared to 
1.9 Mq pc~ 2 Gyr -1 inferred from iKrouna I l|1995l) and 
3 — 5 Mq pc~ 2 Gyr -1 summarized in IR anal 1(1991). Inte- 
grating under the curve, the pressure-based prescription 
gives a total star formation rate for the Milky Way of 
M* = 1.1 M Q yr" 1 compared to 2.0 for the K98 for- 
mulation. Integrating the observed data over the same 
range gives M* = 1.5 Mq yr _1 . 

NGC 598 (M33) - We perform a similar analysis on 
M33 usin g the star fo r matio n rate derived from infrared 
data by IHever et all l)2004[K Again, M33 is entirely 
in the low pressure regime. The integrated value of 
star formation in the galaxy is 0.7 Mq yr _1 in good 
agreement with the v alues found using other methods 
ijHippelein et a l. 2003). In Figure EJ we compare the ob- 
served profile of star formation to that predicted by the 
pressure-based formalism and that of K98. Once again, 
both the pressure-based formalism and that of K98 un- 
derpredi ct the amount of st ar formation. As previously 
noted bv lHever et alJ l|2(304), the K98 prescription does a 
poor job of predicting the slope of the star formation rate 
in M33. The pressure-based prediction agrees with the 
slope of the observational data, only the values are low 
by a factor of seven. Yet Tabled shows that the pressure 
formulation predicts the molecular gas fraction reason- 
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of K98 (dashed) and the m easured values, 
ans IGuesten fc Mezgerl IT982T) : Circles: Pul- 



FlG. 5. — Star formation rate as a function of galactocentric 
radius (R ga i) for the Milky Way. The figure compares the results 
for the pressure-based star formation prescription in the present 
work (solid) to that 
Tria ngles: H II region s _ 

sars ILvne et alj|1985l) : Squares: SNH IGuibert et alj|1978l) . The 

pressure formulation predicts not only the shape of the curve for 
the observed star formation rate, but for most of the points, the 
amplitude as well. 



100.0 



10.0 




Fig. 6. — Star formation rate as a function of galactocentric 
radius (R ga i) for NGC 598 (M33). The figure compares the results 
for the pressure-based star formation prescription in the present 
work (solid) to th at of K98 (dashed) and the measured values (gray, 
IHever et al.l2004ft . The pressure formulation predicts the shape of 
the observed curve, though the results are offset by a constant 
factor of 7, implying the star formation efficiency of molecular gas 
in M33 is dramatically higher than in comparable systems. 



ably well in M33. The discrepancy must arise in the 
efficiency of star formation in the dense gas, e, the dense 
gas fraction, Eh2(HCN)/Eh2(CO) or the calibration of 
either the CO or the Ha data. Future observational work 
will be required to determine why the molecular gas in 
M33 appears to be so efficient at forming stars. 

Molecule-Rich Galaxies - Finally, we examine the pre- 
dictions of the pressure-based star formation pres c riptio n 
in the molecule- rich galaxies of IWong fc Blitzl $2002) . 
In these galaxies, the measured Ejj2(CO) is everywhere 
greater than Ehi, so we use Equation 1251 to determine 
the SFR. We compare the K98 prescription and pressure- 
based results for six galaxies in Figure[7| We omit M101 
(NGC 5457) from our analysis because of the narrow 
range of R ga i over which all relevant quantities are mea- 
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TABLE 3 

Results 



Galaxy M H 2 M+, Observed M*, Pressure Af*, K98 

(10 9 Af o ) (A/^yr- 1 ) (M Q yr" 1 ) (M Q yr" 1 ) 



NGC 4321 6.5 1.8 10.4 6.6 

NGC 4414 4.4 3.0 6.2 4.9 

NGC 4501 1.6 3.2 1.4 1.2 

NGC 4736 0.3 0.3 0.2 0.4 

NGC 5033 4.2 1.9 5.2 3.8 

NGC 5055 1.6 1.2 1.5 1.7 



portant. Here, we argue that the state of the gas is of 
central importance, and can be determined by Equation 
IT51 from the total amount of gas and knowledge of the 
gas pressure. Thus, high column density, low-pressure 
gas should have relatively few molecules and therefore 
a low star formation rate. High column density, high- 
pressure gas should have about the same star formation 
rate in both formulations, because there will be nearly 
complete conversion of atomic to molecular gas. 



5. SUMMARY AND CONCLUSIONS 

• We show that the ratio of atomic to molecular gas 
in galaxies is determined by hydrostatic pressure 
and that the relationship between the two is nearly 
linear, over three orders of magnitude in pressure. 
The rms deviation of individual points from the 
mean relation is a 0.26 dex, or 80%. The scatter of 
the binned data for all of the galaxies is 0.15 dex, 
only 38%. 

• In addition to very similar slopes, the galaxies have 
only a small range in zero points, except for three 
galaxies that are interacting with their environ- 
ments. 

• The dispersion of measurements within individual 
galaxies varies widely from galaxy to galaxy, with 
some galaxies showing extraordinarily small disper- 
sions. 

• Although some of the correlation is driven by the 
appearance of Sh2 on both axes at values of high 
Rmoi, measurement of the power law index of the 
pressure is nevertheless reliable. 

• Although the galaxies in our sample are missing 
SO, E, strong starburst and Seyferts, we probe a 
large range in mass, metallicity and column density, 
suggesting that the P ex t — R m oi relation has wide 
applicability. 

• Even though star formation is associated with 
dense gas traced by CS and HCN, we show that 
CO is as good a tracer of star formation in galaxies 
as HCN, even in starbursts. This is due to the 
tight, non- linear relationship between Ssfr and 
Sh2(CO) in galaxies. 

• We derive a pressure-based prescription for star for- 
mation in galaxies that is based on the P — R mo i 
relation derived in this paper. It differs from K98 
for the outer regions of large spirals, dwarf galaxies, 
damped Lyman alpha systems and other molecule- 
poor galaxies. For molecule-rich galaxies, the pre- 
scription has a form similar to that of K98. 
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sured. The observe d star formation data are those of 
iWong fc Blitl l)2002|) for the case where A v oc N(B). For 
these galaxies, the pressure-based star formation predic- 
tion agrees well with that of K98 since the galaxies are 
predominantly molecular. The two predictions agree rea- 
son ably well with o bserv ed data except for NGC 4321. 
The lWong fc Blit3 l|2002fl observations are insufficiently 
sensitive to go into the low pressure regime where the 
pressure-based formulation may do better at reproduc- 
ing the observed results. Indeed, it is when R mo i < 2 
that the relations diverge, which will occur for the outer 
regions of large spirals, H I-rich dwarfs such as dwarf ir- 
regulars, and damped Ly-a systems. For these objects, 
our predicted star format i on rat es are significantly lower 
than those of iKennicuttl l)1998f) . Since dwarf irregulars 
often show episodic star formation, our prescription as- 
sumes that enough dlrrs are observed to be representa- 
tive of the population as a whole. 

We also compare the observed global star formation 
rates in these Table The global star formation rates 
given in the table are determined by integrating under 
each of the curves. On the basis of the data in this 
table, there is no significant difference in using K98 or 
our Eauation l25l to predict star formation rates for these 
large, molecule-rich galaxies, though the pressure-based 
prescription tends to give somewhat higher values. A 
difference of about 15% is expected as a result of a some- 
what higher value of the coefficient in Equation El com- 
pared to the K98 value. 

galaxies 

4.2.4. Comparison with Simulations 

A fundamental difference between the two relations 
comes when one wishes to compare the predictions of 
the star formation rate from numerical simulations with 
that of the star formation prescriptions. In the K98 for- 
mulations, since the star formation rate depends only on 
the total amount of gas, the state of the gas is unim- 
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